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Abstract 

An  efficient  solution  to  the  3-D  scattering  and  absorption  problems  of  chirality  and  biisotropy  is 
developed  using  various  types  of  expansion  functions.  The  necessary  conditions  to  obtain  the  most 
optimal  solution  are  outlined.  The  program  complex  to  compute  and  visualize  scattering, 
absorption,  polarization,  energetic,  and  directive  properties  of  biisotropic  and  chiral  objects  is 
described  and  applied  to  solve  some  applied  scattering  and  antenna  problems. 

1.  Introduction 

The  problems  of  microwave  behavior  of  the  objects  of  complicated  shape  and  complex  material 
filling,  such  as  chiral  and  biisotropic,  are  to  be  important  both  from  theoretical  and  applied  viewpoints 
[1-6].  In  this  regard,  much  effort  has  been  devoted  in  last  few  years  to  develop  and  justify  new 
techniques  and  codes  to  solve  2D  and  3D  problems  on  chirality  and  biisotropy.  However,  most  of  the 
results  declared  for  the  time  being,  is  restricted  to  the  geometry  of  objects,  the  range  of  material 
parameters,  the  type  of  primary  excitation  and  the  accuracy  attained  and  thus  no  longer  satisfy  the 
present-day  practical  needs. 

In  this  contribution,  based  on  the  Method  of  Auxiliary  Sources  (MAS)  [4-6],  an  efficient 
solution  to  the  3-D  scattering  and  absorption  problems  of  chirality  and  biisotropy  is  developed  using 
various  types  of  expansion  functions.  Farther,  the  necessary  conditions  to  obtain  the  most  optimal 
solution  are  outlined.  The  program  complex  to  compute  and  visualize  scattering,  absorption, 
polarization,  energetic,  and  directive  properties  of  biisotropic  and  chiral  objects  is  described. 
Application  of  the  created  computer  code  to  solve  the  practical  problems  of  interest  is  illustrated. 

2.  Brief  Description  of  the  Method 

The  problem  discussed  here  is  to  find  the  electromagnetic  response  of  an  arbitrary  homogeneous 
biisotropic  (including,  chiral)  object  to  be  excited  from  outside  or  within  by  the  given  field  of  primary 
electromagnetic  sources.  From  the  mathematical  point  of  view,  this  problem  can  be  obviously  reduced 
to  solving  wave  equations 

3 e’‘U(r)  = 0 (1) 

for  unknown  potential  function  of  scattered  field  U(r),  with  providing  the  radiation  conditions  at 
infinity  and  satisfying  the  following  boundary  conditions  on  the  surface  5 enclosing  the  object 

W $(r)+U0(r) Jjr  = r»  = 0,  M(r*  )e  S (2) 

Hereinafter,  superscripts  e and  i concern  to  the  exterior  and  interior  domains,  respectively,  U0(r)  is  a 
given  function  of  incident  field,  and  W is  an  operator  of  boundary  conditions  ensuring  the  fields 
conjugation  on  the  boundary  surface  S. 

Various  techniques  to  solve  the  boundary  problem  (l)-(2)  differ  actually  by  the  method  to  construct 
the  solution  to  the  equation  (1)  and  to  determine  the  unknown  coefficients  arising  from  satisfying  the 
boundary  conditions  (2).  The  MAS,  in  contrast  to  most  other  methods,  constructs  the  solutions  to  (1) 
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as  non-orthogonal  expansions  in  terms  of  fundamental  solutions  of  (1)  with  singularities  out  of  the 
described  domain  [6] 


Nl’e 

U(f)=^a^n(r,r^X  ? e Die  (3) 

n= 1 

r ■ i Ni,e 

The  sets  of  points  ]n=]  to  be  distributed  at  the  auxiliary  surfaces 

Sl'e  (Fig.  1)  can  be  interpreted  as  the  centers  of  auxiliary  sources 
associated  with  fundamental  solutions  Un(r,r„e ).  The  amplitudes  af 
of  these  sources  should  be  determined  from  the  boundary  conditions  (2) 
by  any  numerical  procedure,  in  particular,  by  the  collocation  method. 


3.  Application  to  the  Problems  of  Chirality  and  Biisotropy 


Application  of  the  described  scheme  to  the  problems  of  chirality  and  biisotropy  is  of  the  simplest  one 
if  using  constitutive  relations  in  Post  notations 

D = eE  + iaB,  H = $£  + p.-1  B , (4) 

where  £ and  |1  are  the  medium  permittivity  and  permeability,  respectively,  a = £-i\|t  and 
P = £ + i\|/  are  magnetoelectric  admittances,  £,  is  chirality  admittance  and  \|/is  nonreciprocity 
susceptance. 

To  determine  now  fundamental  solutions  of  wave  equation  (1)  for  constitutive  relations  (4),  it  is 
necessary  first  choose  the  type  of  the  potential  function  U . As  such  a function,  any  function  can  be 
employed  which  identically  define  the  vectors  of  electromagnetic  field,  e.g.,  electric  and  magnetic 

vector  potentials  Ae'm , Debye  scalar  potentials  e’mIl,  Hertz  vectors  Ze,m , spinor  dyad  of  Hertz 
potential  Z , spinor  dyad  of  electromagnetic  field  F or  any  field  vector  {e,  H,  D,  b\.  It  is 
convenient,  however,  to  choose  such  a potential  function,  which  leads  to  the  simpler  form  of  wave 
operator  3 and  the  relations  to  the  field  vectors. 

Thus,  when  choosing  E or  H as  potential  function  U , wave  operator  3 is  of  the  form 

3=(V2+fc2)7  + co}i(a  + p)Vx7,  k = (nje\i,  (5) 

and  fundamental  solutions  of  wave  equation  (1)  are  found  to  be  as  follows 

Vfr,  n = [yrGr(\r  -?])- yec£(\r  - r]  JsRT(r')  , (6) 

where  yrJ  = V x 7 ± ± kr  l I are  vector  differential  operators  for  right-hand  (r)  and  left-hand 


„ r — 

(£)  wavefields,  91  =--l  -pg-  is  normalized  coefficient,  7 is  a unit  matrix, 

con 


Gr’e  =- 


ik 


>.« 


4ujr  - 

are  scalar  Green’s  functions  for  outgoing  waves  with  wavenumbers 

kr,e  =k{[J  + t|2(a  + P)2  /471/2  ±ri(a  + pJ/2;,  t\  = J\Ue, 


and  wave  impedances 

r\rc'£  =T\/{l  + r\2(a  + $)2 /4]112  ±T\(a-$)/2}, 
%(r')  is  a unit  vector  tangential  to  the  auxiliary  surface  at  the  point  M'(r') . 


(7) 

(8) 


(9) 
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If  introducing  as  potential  function  the  spinor  dyad  of  electromagnetic  field  U = F , wave 
operator  3 is  of  the  form 

$ = Vx7-ifc7  (10) 

and  the  fundamental  solutions  of  (1)  are  obtained  to  be  as  follows 

'F(r-r')  = yG(r-r')9lf  (11) 


where  k 


fkr 

0 


-k, 


, Y = 


'r 


0 


0 -f 


f 


and  G = 


G ' 
0 


0 

G ' 


are  the  matrices  of  wave-numbers, 


vector  differential  operators  and  Green's  functions  for  right-  and  left-hand  wavefields,  respectively. 
The  field  vectors  E and  H can  be  found  over  the  potential  function  via  relations 

E = XF,H=ZF,  (12) 

with  rows  of  parameters 

X = (ir\rc  ),  C=(U)  (13) 

The  knowledge  of  fundamental  solutions  of  wave  equation  (1)  for  biisotropic  (including,  chiral) 
medium  formally  completes  the  construction  of  solution  (3)  to  the  stated  boundary  problem  (l)-(2). 


4.  Numerical  Results 


To  numerically  calculate  electromagnetic  characteristics  of  chiral  or  biisotropic  objects  by  the  present 
method,  it  is  necessary  to  determine  the  coefficients  a‘r;e  of  expansions  (3)  from  the  boundary 
conditions  (2).  To  provide  the  quick  convergence  and  minimal  discrepancy  of  the  solution  to  be 
sought,  the  proper  choice  of  the  auxiliary  parameters,  such  as  the  shape  and  dimensions  of  auxiliary 
surface,  distribution  of  the  points  on  the  boundary  and  auxiliary  surfaces,  etc.,  is  required.  It  is 
especially  important  to  properly  account  for  the  character  and  situation  of  the  singularities  of 
continuously  extended  scattered  field  across  the  boundary  of  the  domain  to  be  described. 

All  the  aspects  above  have  been  considered  to  produce  the  program  complex  for  calculating  and 
visualizing  scattering,  absorption,  polarization,  energetic,  and  directive  characteristics  of  3-D  chiral 
and  biisotropic  objects.  In  the  code  performed,  both  isolated  and  integrated  auxiliary  sources  have 
been  utilized  to  obtain  the  optimized  and  effective  solutions.  We  illustrate  here  some  possibilities  of 
the  created  code  to  calculate  and  visualize  radiation,  propagation  and  scattering  problems  of  interest. 

Fig.  2 and  3 compare  the  radiation  processes  produced  by  a linear  current  source  placed  in  a 
back  focus  of  biisotropic  ellipsoid  of  revolution  along  (Fig.  2)  and  perpendicular  (Fig.  3)  to  the  axis  of 
revolution.  The  geometry  of  ellipsoid  is  determined  by  its  semi-axises:  a = 0.5,  b = 0.4,  and  medium 
parameters  are:  £r  =2.00006,  (Xr  =2.02502,  ar|0  =0.1885-1-/0.754,  [3r|0  =0.1 885- i 0.754 

(r|0  =120n).  We  are  interested  in  distribution  of  energy  density  of  right-hand  (a)  and  left-hand  ( b ) 
wavefields  outside  the  ellipsoid  for  wavenumber  k = 32.  The  parameters  of  biisotropic  medium  are 
chosen  so,  that  mainly  left-hand  polarization  is  focused  along  the  axis  of  ellipsoid.  The  comparison  of 
Fig.  2 and  3 also  shows,  that  the  transverse  current  source  manifests  more  focusing  properties. 

Fig.  4 shows  the  3-D  radiation  pattern  of  ellipsoid  of  revolution  of  Fig.  2 and  3 with  longitudinal 
(a)  and  transverse  ( b ) excitation.  From  comparison  of  Fig.  4-a)  and  4-b),  it  is  obvious,  that  transverse 
current  excitation  leads  to  better  directive  pattern  with  weaker  lobes.  It  is  also  clear,  that  this  pattern  is 
not  omnidirectional,  and  thus  one  of  polarization  (left-hand)  is  dominant.  Fig.  5 calculated  for 
situation  of  Fig.  4-a),  but  with  k = 400,  shows  the  increase  of  directivity  and  decrease  of  radiation 
lobes  with  growing  wavenumber  (frequency). 

It  should  be  noted,  that  Fig.  2-5  have  been  calculated  with  rather  high  accuracy  (about  0.5%) 
and  only  thus  allow  studying  in  detail  electromagnetic  properties  of  scattering  biisotropic  objects. 
Thus,  the  proper  use  of  the  method  described  allows  one  to  obtain  numerical  results  with  predesigned 
accuracy  and  to  study  the  practical  problems  of  interest. 
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